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Abstract. - The metal-insulator transition (MIT) of VeOn is studied by means of electronic 
structure calculations using the augmented spherical wave method. The calculations are based 
on density functional theory and the local density approximation. Changes of the electronic 
structure at the MIT are discussed in relation to the structural transformations occuring si- 
multaneously. The analysis is based on a unified point of view of the crystal structures of 
VeOn, VO2, and V2O3. This allows to group the electronic bands into states behaving similar 
to the dioxide or the sesquioxide. While the sesquioxide-like V 3d yz states show rather weak 
changes on entering the low-temperature structure, some of the dioxide-like V 3d x 2_ y 2 states 
display splittings and shifts similar to those known from VO2. The MIT of VeOn arises as 
a combination of changes appearing in both of these compounds. Our results shed new light 
onto the role of particular electronic states for the MIT of V2O3. 



The binary oxides of vanadium are attracting considerable interest for many years. This 
is due to metal-insulator transitions (MIT) displayed by the majority of these compounds, 
in particular their sensitivity to both electronic correlations and electron-phonon coupling. 
The delicate interplay of both mechanisms prevents a straightforward understanding of the 
transitions. As a consequence, the origins of the transitions are still a matter of dispute. This 
holds especially for the prototypical compounds VO2 and V2O3 

The MIT of VO2 at 340 K is accompanied by a structural transformation from the rutilc 
structure to a monoclinic structure. While both electron-lattice interaction and electronic 
correlations have been proposed as the origin of the transition (l],||,||], recent electronic struc- 
ture calculations for VO2 gave strong hints at a Peierls instability of the one-dimensional rfii 
(d x 2_ y 2) band in an embedding background of the remaining V 3d ti g states |^,|(|. This sce- 
nario was supported by studies of the neighbouring compounds M0O2 and NbC>2 Yet, 
since calculations based on the local density approximation (LDA) just miss the insulating 
gap of low-temperature VO2, electronic correlations are also expected to be of some relevance. 

At a temperature of 168 K and ambient pressure, stoichiometric V2O3 experiences a transi- 
tion from a paramagnetic metallic (PM) phase to an antifcrromagnetic insulating (AFI) phase. 
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At the same time, the crystal structure transforms from the corundum to a monoclinic struc- 
ture. In contrast, on doping with only small amounts of Al or Cr, V2O3 enters a paramagnetic 
insulating (PI) phase without change of the crystal symmetry. Yet, recent EXAFS experi- 
ments gave evidence of local crystal structure distortions in the PI phase, which only in the 
monoclinic AFI phase display long-range order The metal-insulator transitions of V2O3 
are generally regarded as classical examples of a Mott- Hubbard transition |HJ . LDA calcula- 
tions performed for all three phases showed only a minor response of the electronic structures 
to the crystal parameter changes. The slight narrowing of the characteristic a\ g bands in 
the insulating phases formed the basis for a successful description of the PM-PI transition 
in a recent study, where LDA calculations were combined with dynamical mean-field theory 
(DMFT) (jnj. This approach, called LDA+DMFT, is especially suited for a realistic modelling 
of strongly interacting electron systems. The strong influence of electronic correlations has 
also been demonstrated by LDA+U calculations for the AFI phase jl2| . 

A more comprehensive understanding of both compounds calls for deeper insight into the 
role of electronic correlations and electron-lattice interaction on going from VO2 to V2O3. 
In particular, we are interested in the electronic states which are involved in the MIT. In 
the present paper we address these issues by discussing the results of electronic structure 
calculations for the Magneli phase VgOn. The Magneli phases of vanadium form a homologous 
series V„02 n -i (3 < n < 9) and are particularly suited for studying the differences in crystal 
structures and electronic properties between the end members VO2 (n —> 00) and V2O3 
(n = 2). Being part of a broad investigation of the Magneli phases, our work is the first 
ab initio study of these compounds at all. Relating the different local environments of the 
vanadium atoms in VgOn to the electronic properties, we are able to obtain deeper insight 
into the crossover from the dioxide to the sesquioxide. In particular, we find that the changes 
of the d x 2_ y 2 states along the series account for most of the differences. 

As the general formula V„02 n -i = V2O3 + (n — 2)VC>2 suggests, the crystal structures 
of the Magneli phases are usually viewed as rutile-type slabs of infinite extension, which are 
separated by shear planes with a corundum-like atomic arrangement 

[[!,§, [U|- While in the 

rutile-type regions the characteristic VC>6 octahedra are coupled via edges, the shear planes 
have face-sharing octahedra as in the corundum structure. 

For a deeper understanding of the electronic properties of all Magneli phases we have 
developed a different representation of the crystal structures. It starts out from the regular 
3D network of oxygen octahedra, which, apart from a slightly different buckling, is similar for 
all members including VO2 and V2O3. Differences between the compounds V n 02n-i arise 
from the filling of the octahedra with vanadium atoms. Filled octahedra form chains of length 
n parallel to the pseudo-rutile axis c pru t, followed by n — 1 empty sites. In particular, in 
VO2 these chains have infinite length. For VgOn the situation is sketched in Fig. [[[ which 
shows a projection of an O-V-0 sandwich-like slab cut out of the crystal structure. In this 
projection, oxygen octahedra appear as a regular hexagonal network. Both this network and 
the vanadium chains of length 6 can be identified in the left panel of Fig. [I]. Perpendicular 
to the projection, along a pru t, vanadium and oxygen layers alternate. Two different types of 
vanadium layers are distinguished, which comprise vanadium atoms VI, V3, V5, and V2, V4, 
V6, respectively. Within the vanadium sublattice these layers alternate along a pru t- However, 
vanadium chains in neighbouring vanadium layers are shifted parallel to c pru t, such that atoms 
VI and V2 are found on top of each other. This situation is sketched on the right panel of Fig. 
[l]. The relative shifts of the vanadium chains give rise to different arrangements of the VOg 
octahedra. While octahedra neighbouring along a pmt or 6 prut share faces, coupling between 
metal atoms along all other directions within the layers is via octahdral edges. Whereas the 
atomic arrangements near the ends of the chains are similar to those in V2O3, coordination 
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Fig. 1 - Crystal structure of VsOn. Left: View along a pru t. Large and small circles mark vanadium 
and oxygen atoms, respectively. Right: V-chains V1-V3-V5-V5-V3-V1 (black circles) and V2-V4- 
V6-V6-V4-V2 (gray circles) belonging to adjacent vanadium layers. The table lists measured V-V- 
distances (A) 113] . 



near the chain centers is the same as in the rutile structure. The projection of the latter is 
indicated by a rectangle. Our representation offers the great advantage of allowing reference 
of the single symmetry components of the V 3d orbitals to coordinate systems, which are 
common to the whole class of compounds. Inspired by previous work on the dioxides we will 
use local coordinate systems with the z and x axis parallel to the apical axis of the local 
octahedron and the pseudo-rutile c pru t axis, respectively 

The MIT of VgOn at 170 K is accompanied by a structural transformation, which preserves 
the triclinic space group PI. Major changes concern the strong V4-V6 dimerization as well as 
ferroelectric-like shifts of atoms V3, V4, V5, and V6 away from the centers of the surrounding 
oxygen octahedra. In contrast, the V5-V5 dimerization and the strong ferroelectric-like shifts 
of atoms VI and V2 are present in both phases with only minor changes. Hence, even with 
respect to the details of the local octahedral distortions, the chain centers and ends behave 
similar to their parent structures VO2 and V2O3. 

The LDA calculations were performed using the scalar-relativistic augmented spherical 
wave (ASW) method jl6| with a standard parametrization of the exchange-correlation poten- 
tial |l7|. Crystallographic data given by Canfield @ were used. In order to account for the 
openness of the crystal structures, socallcd empty spheres were included to model the correct 
shape of the crystal potential in large voids. Optimal empty sphere positions and radii of 
all spheres were determined automatically |Q. The basis set comprised V 4s, Ap, 3d and O 
2s, 2p as well as empty sphere states. Brillouin zone sampling was done using an increased 
number of k points ranging from 108 to 2048 points within the irreducible wedge. 

Partial V 3d densities of states (DOS) resulting from calculations for the crystal structures 
of both phases are displayed in Figs. || and ||. The gross features are very similar to those 
known from VO2 and V2O3. Two groups of bands are identified, which result from crystal field 
splitting of the V 3d states due to the oxygen octahedra. Note that all orbitals are referred to 
the above local coordinate systems. While V 3d t 2g bands dominate in the energy region from 
—0.7 to 1.8 eV, the e g states (not included in Figs. || and ||) are found mainly between 2.0 and 
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Fig. 2 - Partial VI and V2 3d DOS (per atom) of high- (top) and low-temperature (bottom) VeOn. 

4.8 eV. Small contributions from the ti g states in this energy region result from distortions of 
the octahedra. Not shown are O 2p states occupying the energy region from —8.5 to — 2.5eV. 
Contributions of states in the energy range where the respective other orbital dominates are 
below 10%, but still indicative of covalent bonding. The electronic states obtained for the 
low-temperature phase miss the insulating gap, what is due to the well known limitations of 
the LDA (and very similar to the situation reported for V2O3). However, these limitations 
are not important for the present work, as it aims primarily at the relation of the relevant 
orbitals to different local environments and their changes at the transition. 

Going into more detail, we turn to the t2 g orbitals of those atoms which are involved in 
dioxide-like displacements, i.e. V4, V5, and V6. Their partial DOS are very similar among 
themselves and to those known from the rutile and monoclinic phase of VO2 |§. In particular, 
note the two peak structure of the high-temperature V4 d x 2_ y 2 DOS, which is due to strong 
cr-type metal- metal bonding parallel to the vanadium chains (c pru t). In the low-temperature 
structure, V4-V6 dimerization causes increased splitting of this DOS into bonding and an- 
tibonding branches. At the same time, both the d yz and d xz partial DOS experience an 
energetical upshift due to increased p-d overlap arising from the ferroelectric displacement of 
the vanadium atom perpendicular to c pru t. As in VO2 energetical separation between both 
types of bands is increased albeit not complete. The partial DOS of atom V5 display these 
features already in the metallic phase. In total, the partial DOS of the atoms V4, V5, and V6 
can be fully understood in terms of the mechanisms known from VO2. 

The situation is different for the atoms VI and V2, which are in a sesquioxide-like envi- 
ronment. In particular, they are involved in metal-metal bonding via shared octahedral faces 
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Fig. 3 - Partial V3 and V4 3d DOS (per atom) of high- (top) and low-temperature (bottom) VeOn. 



across the layers. Furthermore, these atoms do not experience any dimerization. As a conse- 
quence, except for small peaks and shoulders near —0.5 and between 1.0 to 1.3 eV, which are 
reminiscent of the respective peaks of the chain center atoms, the d x 2_ y 2 DOS of atoms VI 
and V2 consist of single rather broad peaks extending from —0.2 to 0.7 eV. In both phases, 
subpeaks observed in this energy region do not have counterparts in the DOS of neighbouring 
atoms indicative of metal-metal bonding. Like the d xz states of these atoms the d x 2_ y 2 states 
may thus be regarded as rather localized. This is in contrast to the d yz states, which show 
bonding-antibonding splitting due to the aforementioned overlap across the layers, i.e. parallel 
to ctprut. Since this direction is parallel to the hexagonal Ch ox axis of the corundum structure, 
the splitting of the d yz states is equivalent to that of the a\ g states of V2O3. The partial DOS 
of atoms VI and V2 thus can be interpreted in full accordance with the sesquioxide. 

Atom V3 is exceptional since it does not fit into the above schemes: Neither does this atom 
experience dimerization along the chains as in VO2, nor does it have neighbouring vanadium 
atoms across the layers as in V2O3. Yet, the partial DOS are rather similar to those of atoms 
VI and V2. The d x 2_ y 2 and d xz partial DOS form rather broad single peaks, the former still 
showing satellite peaks at —0.4 and 1.0/1.4eV reminiscent of the bonding and antibonding 
peaks of atom V5. In contrast, the d yz partial DOS displays a pronounced double peak 
structure indicative of bonding-antibonding splitting. This puzzling situation is resolved by 
taking into account a different type of metal-metal bonding, namely that across octahedral 
faces along the sequences V1-V5-V3 and V2-V6-V4 parallel to & pru t- Hence, this cr-type 
overlap and the resulting bonding-antibonding splitting of the d yz states is present not only 
for atom V3, but for all atoms. It is also observed in VO2 f§. For the d yz states of atoms VI 
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and V2 it adds to the above bonding parallel to a prut ; note that the respective antibonding 
peaks merge in the high-temperature phase. Although the coordination of atom V3 differs 
from that of all other atoms, its electronic properties are thus similar to that of the chain 
end atoms due to er-type metal-metal overlap via octahedral faces within the layers. Hence, 
metal-metal bonding across the layers is not the only source of the splitting of the d yz states. 

In conclusion, while overlap of O 2p and V 3d states places the V 3d t2 g states near 
the Fermi energy, the detailed electronic properties of the latter and the MIT of VeOn are 
strongly influenced by the local metal-metal coordination. For this reason, the partial DOS 
for each atom may be regarded as local quantities. Near the chain centers (atoms V4 to V6) 
dimerization and ferroelectric-like displacements via strong electron-lattice interaction cause 
a splitting of the d x 2_ y 2 states and an energetical upshift of the d yz and d xz states. This is in 
complete analogy to the behaviour known from VO2. 

In contrast, near the chain ends (atoms VI, V2, and V3) the d x 2_ y 2 states experience less 
metal-metal overlap and, hence, become much more localized; their partial DOS resemble those 
of the d xz states in the chain centers. As a consequence, electronic correlations might play 
a greater role. Being subject to metal-metal bonding only within the layers, the d yz partial 
DOS of atom V3 displays bonding-antibonding splitting very similar to that of atoms VI and 
V2. In addition, these atoms experience metal-metal overlap across the layers. We conclude 
that the V-V overlap within the layers is at least as important as the one perpendicular. In 
summary, the MIT of VgOn is interpreted as resulting from electron-lattice interaction and 
electronic correlations in the dioxide- and sesquioxide-like regions of the crystal, respectively. 

Our findings have important implications for the understanding of V 2 3 . According to the 
above results the electronic properties of this material are influenced by (i) the localization of 
the d x 2_ y 2 states, and (ii) the splitting of the d yz states due to metal-metal bonding within 
layers perpendicular to the hexagonal Chex &X1S. Splitting of the ai 9 -like states due to overlap 
parallel to Ch cx might be less important than commonly assumed. Both effects could not be 
analysed in previous studies of V2O3 due to the high symmetry of this material. Yet, they 
are accessible to our present analysis of the sesquioxide-like regions of VgOn. Nevertheless, 
the electronic structure as well as the MIT of VgOn remain a crucial test case for all theories 
aiming at a correct description of both VO2 and V2O3. 

* * * 
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